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We present multi-epoch global VLBI observations at 2 cm, 3.6 cm and 6 cm of the compact radio 
sources in Arp220. We resolve many sources and estimate sizes, expansion velocities and source 
classes. We find most source properties are consistent with them being either radio supernovae 
(SNe) or supernova remnants (SNRs). We extend the luminosity-diameter relation for SNRs 
to very small sources and argue this supports models where shell magnetic fields are internally 
amplified. We also detect one probable SN/SNR transition object candidate and three highly 
variable sources with possible superluminal motion (approximately 4c) of jet-like features near 
rapidly varying almost stationary components. These enigmatic sources, which show similarities 
to the recently discovered superluminal source in M82, might be associated with an AGN or a 
new radio source class (e.g. intermediate-mass black holes or micro-blazars). 
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1. Introduction 

1.1 Arp220 

Arp220 is the result of the collision of two galaxies now in the process of merging and has 
two distinct nuclei separated by approximately 350 pc |JIJ]. It is also the closest ULIRG (Ultra 
Luminous Infra-Red Galaxy) from the Earth, located 77 Mpc away. It has a far infrared luminosity 
of Lfir ~ 1.3 x 10 12 L [g] which is more typical of star-forming galaxies at redshift z = 1 and has 
a similar star-formation density per unit area as redshift z = 6 proto-galaxies. Smith et al 1998 [||] 
made the first detection at 18 cm of compact objects associated with the starburst. These objects 
were first detected at higher frequencies (13 cm, 6 cm and 3.6 cm) by Parra et al 2007 [Q] who 
concluded that they are likely a mixture of SNe and SNRs. 

1.2 Supernovae and supernova remnants 

In order to classify the sources observed in Arp220 it is important to understand the differences 
between SNe and SNRs and their radio emission mechanisms. The radio emission of a SN is the 
result of the interaction of the blast wave from the SN explosion with the ionised circumstellar 
medium (CSM). Relativistic particles in the SN blast wave interacting with a magnetic field produce 
non-thermal synchrotron emission which is absorbed by foreground thermal electrons. The amount 
of absorption decreases with time as the distance between the traveling shock and the edge of the 
CSM wind-blown bubble decreases. The radio light-curves characteristically peak from higher to 
lower frequencies followed by a power law decline at each frequency. 

When the shock wave reaches the edge of the wind-blown bubble it starts interacting with 
the denser interstellar medium (ISM) gas. The SN becomes then a SNR. In a first phase the SNR 
expands at a nearly constant velocity radiating thermal X-rays and synchrotron emission. The light- 
curves increase at all frequencies because of the acceleration of relativistic electrons at the shock 
front created where the SNR meets the ISM. When the mass of ISM swept up by the shock wave 
is the same as the initial mass of the stellar ejecta the remnant enters the Sedov-Taylor phase and 
begins to decelerate. At this point it reaches its maximum radio luminosity. There exist different 



theories describing the radio emitting behavior of the SNR from this point on (see § 2.2). 



2. HSA 2 cm and global VLBI 3.6 cm observations 

2.1 Source sizes and classes 

High-resolution HSA 2 cm and global VLBI 3.6 cm observations of Arp220 performed in 
December 2006 resulted in the detection of 17 sources in the Eastern and Western nuclei; a number 
of these sources are resolved - 13 at 2 cm, 8 at 3.6 cm, and 7 at both wavelengths (Batejat et al 
201 1 in preparation). We find a good correlation between source sizes measured at 2 cm and 3.6 cm 
(Fig [[]). Based on their time variability and spectra we classified each source as a SN or a SNR 
(see § [L2| ). All SNRs detected at high frequencies are resolved with sizes d > 0.27 pc. Only one 
source classified as being a SN (E14) is resolved at both bands with a size d = 0.30 pc. All other 
sources have sizes d < 0.22 pc. Sizes are calculated combining the results of our observations at 
both bands. We conclude that our SNe have upper limits on their expansion velocities of 9000 to 
80000 km s" 1 . 
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Figure 1: Source diameter at 2 cm and 3.6 cm with error bars. Red sources are resolved at both bands, green 
only at 2 cm, blue only at 3.6 cm and white are unresolved. Red sources are classified as SNRs except E14. 



2.2 Luminosity - Diameter relation 

Fig ^| plots the detected Arp220 source sizes together with the SNRs detected in M82 and 
in the LMC. The best fit to our observations is a power law relationship between the luminosity 
and the diameter of the SNRs (L oc D~ 19 ) close to the theoretical relation for SNRs in the Sedov 
phase as derived by Berezhko and Volk [[7|] (L oc D~ 9 / 4 ) which also fits our data within the scatter 
(green line in Fig |2|). This model is based on the assumption of an internal B-field energy being 
a constant small fraction of the post shock pressure (pisM-Vf) and therefore declining with the 
source size as the SNR expansion decelerates (see for the details). The data shown in Fig || 
provides strong support for the existence of the Luminosity-Diameter relation with an exponent 
close to that theoretically predicted (Batejat et al 2011 in preparation). Alternative models, such 
as the one presented by Thompson et al [|I0|], claim a constant internal B-field in the Sedov phase 
being simply proportional to the ISM B-field increased by a factor 3 to 6 and predict a constant 
radio luminosity versus size relation. This is inconsistent with our observed Luminosity-Diameter 
relation. Additionally in conflict to this theory the 18 cm flux of individual SNRs decline by 
~ 4%/yr (Batejat et al 201 1 in preparation) which is consistent with Berezhko and Volk [^] model 
if the diameter of the sources expands by 2%/yr (as expected in Arp220 given their measured sizes 
and expected expansion velocities) confirming the idea of a B-field changing with source size. 

3. Recent global VLBI 6 cm monitoring 

We performed deep global VLBI 6 cm observations at three epochs on June 2008, October 
2008 and February 2009 with the goal of looking for short time variability sources of all kinds, i.e. 
due to type Ib/c supernovae or AGN candidates, as well as producing the deepest 6 cm map to date 
of Arp220 by combining the three epochs with a noise rms of about 8/iJy/beam (Fig [3]). 

3.1 New detections 

Nine sources are detected for the first time; five in the Eastern nucleus and four in the Western 
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Figure 2: Illustration of the 3.6 cm luminosities and sizes of observed SNe and SNRs. The SNe tracks 
for SN1986J (Type Iln), SN1979C (Type EL) and SN1980K (Type IIL) were produced using the light 
curve fits given by Weiler et al [|[ combined with the deceleration parameter from Bietenholz et al |Q| for 
SN1986J and assuming free expansion at 10 4 km s _1 for both SN1980K and SN1979C. The square markers 
along each track indicate time evolution and are 1 year apart. Red points = sources identified as SNe, 
blue points = sources identified as SNRs, green point = possible transition object, cyan stars = unclassified 
sources. The green line is the theoretical relation of Berezhko and Volk [0] of slope -9/4. The blue line is 
the best fit considering just those sources identified as SNRs. The data for the SNRs in M82 are taken from 
Huang et al [|] and the data for the SNRs in the LMC from Mills et al [§]. 



nucleus of Arp220. One of these (W61) seems to be a new rising source that was not present at the 
same flux density at earlier 6 cm epochs. 

3.2 Transition object candidate 

Most detected objects are classified as SNe and SNRs. However this campaign of observations 
confirm the unusual characteristics of a few sources. One of them, referred as W12, has light 
curves which increase at all frequencies (3.6 cm, 6 cm and 18 cm) between 2004 and 2009. For 
SNRs we expect declining light-curves at all frequencies and for SNe light curves that rise at low 



frequencies and decline at high frequencies (see § |1.2| ). The observed behavior of W12 makes it a 
good candidate for being a SN/SNR transition object in which the shock front has just reached the 
ISM and electrons are being accelerated giving a rising flux density at all frequencies. 

3.3 Variable sources 

We also have detected three rapidly variable sources. Those three sources, W7, W26 and W29, 
all localised in the Western nucleus of Arp220, are rapidly varying in flux at 6 cm and 1 8 cm, but 
also in position and shape with secondary features appearing and disappearing from epoch to epoch 
(Fig f|). If emitted half way between epochs the speed estimates of the blob features are superlu- 
minal with speeds from Ac to 12c, similar to the recently discovered source in M82 by Muxlow et 
al [O]. The sources W7 and W26 are long-lived sources observed for more than 15 years. W29 is 
greater than 8 years but it is likely much older because it would not have been detected at the same 
flux density in earlier less sensitive epochs. All three sources show behaviours inconsistent with 
conventional SNe or SNRs. It is possible that at least one of those sources is associated with a 
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Figure 3: Greyscale shows the combined 5 GHz global array GC031A,B,C natural weighted image of the 
Western nucleus. The rms noise is 8/xJy/beam, the image size is 512 x 256 mas. Labeled in white are the 
transition object candidate W12 and the highly variable sources W7, W26 and W29. The green crosses give 
the 1.3 mm and 2.3 mm IRAM positions of the hot dust "AGN" feature from Downes and Eckart JIT[], the 
green circle is the position from the 0.8 mm observation of Sakamoto et al jl2|]. Errors on these positions are 
estimated to be of order 50-100 mas. The black ellipse gives the orientation and size of the hot dust feature 
as fitted by Downes and Eckart [ pT| ] centred at the centroid of the three mm interferometer positions. 



control 




Figure 4: Variability of W7, W26 and W29 over time at 5 GHz. The columns show in order the varying 
sources W7, W26 and W29 and finally as a control the non-varying resolved source W18 at epochs A to C 
(from top to bottom). The contours are at 4, 6, 8, ... times the rms noise for the variable sources and at 4, 
12, 20, ... times the rms noise for W18. The rms noise is approximately 15 /iJy/beam at all 3 epochs. The 
beam is plotted for each source on the last row. 



supermassive black hole (SMBH). Downes and Eckart [11] presented dynamical evidence for such 
a SMBH in Arp220's Western nucleus by analysing the position- velocity diagram from CO(2-l) 
observations. Downes and Eckart [JTTJ] and Sakamoto et al [|12|] have analysed the mm continuum 
emission from hot dust and have argued that it may be too compact to be from a starburst and could 
be black hole powered instead. When comparing the position and size of the hot-dust feature from 



Downes and Eckart Jl 1[] and Sakamoto et al [12] with the positions of our three variable objects 
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(see Fig ||) we see that W26 and W29 are located close to the centre of the hot dust emission mak- 
ing them possible AGN candidates. It is however difficult to come to a final conclusion on the 
true nature of those three variable sources. Are those objects accreting supermassive black holes? 
Could they be highly beamed stellar-mass black holes ('micro-blazars') or even radio-loud versions 



of accreting intermediate-mass black holes ||14p. Scheduled HSA observations in the spring 2011 
looking for variability on timescales shorter than 30 days will hopefully give us more clues as to 
the nature of these objects. 

4. Conclusion 

The two observation campaigns presented in this work allow us to conclude that the compact 
radio sources detected in Arp220 are a mixture of SNe and SNRs. We also confirm the theoret- 
ical Berezhko and Volk [Q] Luminosity-Diameter relation for SNRs and extend it to very young 
sources. We have possibly detected a SN/SNR transition object and we have detected at least three 
highly variable sources whose structures vary on timescales shorter than 4 months with possible 
superluminal motion. 
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